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\The  problem  of  simultaneous  heat  and  mass  transfer  in  a  homogeneous  snow  layer,  with  one  side  kept  at  its  initial 
temperature  and  the  other  side  with  a  step  temperature  increase,  was  solved  for  the  case  of  constant  through-flow  condi¬ 
tions.  An  experimentally  determined  effective  thermal  conductivity  function,  i.e.  /Cfc  =  0.0014  +  0.58  G  (where  G  is  dry 
mss  flow  rate  of  air  in  g/cm2s),  was  employed  in  the  solution.  The  computed  nondimensional  temperature  distribution 
agreed  quite  well  with  experimental  data  taken  under  pseudo-steady  state  conditions  with  the  exception  of  the  temper- . 
fature  for  the  lowest  flow  rate  used  in  the  experiment  The  pronounced  nonlinearity  of  the  temperature  distribution  was 
found  to  be  a  strong  function  of  the  flow  rate.  For  sinusoidal  variation  of  atmospheric  pressure,  the  responding  flow  in 
the  snow  medium  was  also  found  to  be  sinusoidal.  In  conjunction  with  the  diurnal  temperature  change,  this  variation  / 


FORM 
t  JAM  7» 


EDITION  OF  I  MOV  IS  IS  I 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  fWNwi  Dott  Rnttrtd) 


SECURITY  CLASSIFICATION  OF  THIS  PAOSfWlln  DM*  MnfndQ 


20.  Abstract  (Corn’d) 

facilitated  the  process  of  repeated  sublimation  and  condensation  in  alternate  directions  and  thereby  produced  a  surface 
layer  of  approximately  constant  snow  density. 


Unclassified 

security  classification  of  this  RAOCprfiMi  DM*  But* r*« 


PREFACE 


This  report  was  prepared  by  Dr.  Yin-Chao  Yen,  Chief,  Physical  Sciences  Branch,  Research 
Division,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory.  Funding  for  this  re¬ 
search  was  provided  by  DA  Project  4A1 61 102AT24,  Research  in  Snow,  Ice  and  Frozen  Ground, 
Task  A,  Properties  of  Cold  Regions  Materials,  Work  Unit  005,  Thermophysics  and  Mechanics  of 
Cold  Regions  Materials. 

The  manuscript  of  this  report  was  technically  reviewed  by  Dr.  Samuel  Colbeck  and  Dr. 
Yoshisuke  Nakano  of  CRREL. 


iii 


fQ 


NOMENCLATURE 


heat  capacity 
vapor  diffusivity 
k  permeability 
K  thermal  conductivity 
L  latent  heat 

t  layer  depth 

M  molecular  weight 

p  partial  pressure  of  water  vapor,  vapor  pressure, 
air  pressure 

P  defined  asp2  -p% 

T  temperature 
t  time 

v  defined  by  (7"  -  T0) 
w  defined  by  eq  1 2 

ae  defined  by  eq  8 

P  defined  by  eq  9 


y  defined  by  eq  20 

r)  dimensionless  distance  (y/i);  dimensionless 
number  defined  as  (w/2a)^ 
n  system  pressure;  a  constant  =  3.1 41 6 
p  density 

0  dimensionless  temperature  (T-T0)l(Tt-TQ)-, 
also  as  defined  in  eq  26 
<7  defined  as  kp0/e 
u  period 

Subscripts 
a  air 

e  effective;  porosity 
fi  at  y  =  8 
s  snow;  sublimation 
w  water 
0  aty^O;  initial 


ON  THE  TEMPERATURE  DISTRIBUTION 
IN  AN  AIR-VENTILATED  SNOW  LAYER 


Yin-Chao  Yen 


INTRODUCTION 

In  nature,  the  surface  layers  of  a  snow  cover  are  affected  by  the  constantly  changing  temperature 
of  the  ambient  air.  During  the  winter  months,  the  bottom  layers  have  a  significantly  higher  and 
more  constant  temperature  than  the  atmosphere.  The  temperature  difference  between  the  upper 
and  lower  layers  gives  rise  to  various  processes  within  the  snow  cover  that  are  caused,  in  part,  by 
the  different  vapor  pressures  between  the  snow  grains.  Kondrateva  [1945]  reported  on  the  thermal 
conductivity  of  artificially  compressed  snow.  She  indicated  qualitatively  the  extent  of  the  vapor¬ 
ization-sublimation  processes  occurring  in  the  snow  sample  by  comparing  the  snow  density  of  each 
layer  determined  at  the  end  of  an  experiment  with  those  at  the  beginning.  Yosida  (1950)  experi¬ 
mentally  determined  the  macroscopic  diffusion  coefficient  D  (or  the  effective  diffusion  coefficient 
Dt)  of  water  vapor  through  snow  by  measuring  the  weight  change  of  a  snow  sample  under  a  con¬ 
stant  temperature  differential.  (This  coefficient  includes  all  the  possible  mass  transfer  mechanisms 
occurring  in  a  snow  medium,  i.e.  migration  of  water  molecules  along  the  snow  grains,  sublimation 
and  condensation,  and  diffusion  through  the  void  space  between  the  snow  grains.)  He  concluded 
1 )  that  the  diffusion  coefficient  of  water  vapor  through  snow  is  independent  of  snow  density,  2) 
that  it  is  four  to  five  times  larger  than  the  diffusion  coefficient  D0  of  water  vapor  through  air 
(0.22  cm2/s  at  0°C  and  atmospheric  pressure)  and  3)  that  gravity  has  no  effect  on  the  value  of  the 
diffusion  coefficient. 

In  these  two  investigations,  the  snow  samples  were  subjected  to  a  constant  temperature  gradient 
with  stagnant  fluid  (air  in  this  case)  in  the  void  space  between  the  snow  grains.  In  nature,  however, 
a  certain  degree  of  convective  transfer  of  heat  and  mass  occurs  near  the  surface  layer  of  snow  cover 
owing  to  wind  currents.  Yen  (1962, 1963)  studied  the  effect  of  air  flow  on  the  effective  thermal 
conductivity  Kt  and  vapor  diffusivity  Dt  and  found  the  following  relations: 

Ke  =  0.001 4  +  0.589  G  (1) 

and 

Dt  =  95.38  (G  +  0.456x  lO'4)5*  (2) 

for  snow  densities  from  0.376  to  0.472  g/cm3  and  a  mass  flow  rate  of  dry  air  G  from  2x  1 0"4  to 
40x  1 0*4  g/cm2  s.  Adivarahan  (1 961 )  extensively  investigated  the  effective  thermal  conductivity 
of  porous  rocks  under  the  influence  of  fluid  flow  but  without  the  accompanying  effect  of  mass 
transfer. 

The  purpose  of  this  study  is  to  demonstrate  analytically  the  resultant  effect  of  an  imposed 
temperature  differential  and  a  constant  through-flow  of  air  on  the  temperature  distribution  of  an 
initially  homogeneous  snow  layer.  The  results  will  be  compared  with  reported  data  on  flow  through 
consolidated  and  unconsolidated  porous  media  with  or  without  mass  transport  To  consider  the 


r 


more  practical  case  of  variable  flow  situations,  a  sinusoidal  atmospheric  pressure  variation  on  the 
snow  surface  is  considered  and  the  characteristic  nature  of  the  penetrating  wave  into  the  medium 
is  derived.  The  combined  effects  of  diurnal  temperature  change  and  propagating  wave  character¬ 
istics  on  the  processes  of  heat  and  mass  transfer  on  snow  mass  redistribution  are  assessed. 


FORMULATION  OF  THE  PROBLEM 


Constant  flow 

If  we  consider  a  homogeneous  snow  layer  with  constant  countercurrent  or  cocurrent  air  flow 
rates,  the  energy  balance  equation  for  a  differential  layer  thickness  dy  can  be  written  as 


d2T 
dy 2 


+  G  C  —  + 

'  °  Cpa  dy  ~ 


GM  w 


dT 

a  t 


where 
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w 
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pa 


y  =  distance  measured  from  the  direction 
of  heat  flow 

Kt  -  effective  thermal  conductivity 
p  =  partial  pressure  of  the  water  vapor 
rr  =  system  pressure 

and  Ma  =  the  molecular  weight  of  water  vapor 
and  dry  air 

and  Cps  =  the  heat  capacities  of  air  and  snow 
ps  =  snow  density. 


(3) 


This  equation  is  derived  from  the  assumptions  that  1)  the  flow  is  unidirectional  at  a  constant  rate, 
2)  the  density  of  the  air  is  constant,  and  3)  the  temperatures  of  the  air  and  of  the  snow  are  the 
same  at  the  same  point  in  the  medium  and  time.  Based  on  the  work  of  Bader  et  al.  (1939),  the 
partial  pressure  of  water  vapor  in  air  can  reasonably  be  replaced  by  the  saturation  vapor  pressure 
of  snowp,  which,  as  reported  by  Yosida  (1950),  can  be  represented  by 

pf  =  6.1  exp  (0.0857  T)  (4) 

where  pi  is  in  millibars  and  T  is  in  degrees  Celsius.  For  a  temperature  range  of  -1 7°  <  T  <  -7°C, 
eq  4  can  well  be  approximated  (with  an  average  error  of  ±  4.0%)  by 


p  =4.636  +  0.195  T. 


(5) 


Substituting  eq  1  and  5  into  3  and  considering  the  case  that  7r  >>  pv  eq  3  becomes 
(0.0014  +  0.586}  G(c„  +0.1 95  %£)  f  |f 


(6) 


where  the  +  sign  indicates  that  air  flow  is  in  the  opposite  direction  of  the  heat  flow,  and  Li  is  the 
latent  heat  of  sublimation.  The  term  0.195  represents  the  contribution  associated 

with  net  water  vapor  sublimation.  The  initial  and  boundary  conditions  to  be  considered  are 

T(y,0)=T0 

mo  =  r0 
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2 


(7) 


mo  =  Tf 

By  using  the  transformation  of  v(y,  t)  =  T  ~  T0  and  defining 


_  0.0014  + 0.58  G 
"  (CpsP5) 


and 


±  G 


(c»*ol95%^ 


Cp$P 


eq  6  becomes 


a !e  ^+/3^  = 

3y2  3y  3? 


with  corresponding  initial  and  boundary  conditions 
Wy,0)  =  0 
WO./)  =  0 
W« ./)  =  Tt-T0. 

By  using  the  further  transformation  of 
/»2y  1 


v'  =  vv  exp 


i2ae  4<i( 


eq  1 0  and  1 1  become 
32tv  _  3w 


3y2  3r 


and 


wfy.O)  =  0 
M0.r)  =  0 


MM  =  (rs-r0)exp(|-+^). 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


The  solution  of  eq  13  with  the  initial  and  boundary  conditions  of  eq  14  can  be  obtained  through 

_ _ _  r\..i _ _ i  /■  irii  _ i  i  rtr/-\  _ •  »  _• _ _ 


the  use  of  Duhamal’s  formula  (Hildebrand  1956)  and  is  given  as 


.  nny 

-  sin— 

(32E2+4ae27r2n2  ® 
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n- 1  n 


.  nny 
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(15) 
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The  last  term  accounts  for  the  transient  state  behavior.  But  as  time  goes  by  it  becomes  less  signifi¬ 
cant,  and  the  temperature  distribution  for  the  steady  state  is 


(-1  )"02«2  .  nity 

—  ~  -  -  — .  5| p  -  * 

/?(/?2C2  +  4<Jf2ff2(72)  ® 


(16) 


Equation  16  is  hard  to  evaluate  because  of  the  slow  convergence  of  the  series.  However,  in  the 
case  where  02£2  is  much  smaller  than  4<ar  7T2,  the  following  approximation  can  then  be  made: 


A  positive  value  of  0  indicates  that  the  direction  of  air  flow  is  opposite  that  of  heat  flow  (or  coun¬ 
tercurrent  flow),  and  a  negative  value  of  0  shows  the  same  direction  of  heat  and  air  flow  (cocurrent 
flow). 

For  the  steady-state  case,  eq  6  is  a  linear  second-order  differential  equation  with  constant  co¬ 
efficients,  and  its  solution  with  its  initial  and  boundary  conditions  shown  in  eq  7  can  be  written  as 


T_h__  1  -exp  (-7 y) 
Ts-T0  l-exp(-yS) 


(19) 


where  7  is  defined  as 

y  -+G  Cpa  +  °-19s 

'  0.0014  +  0.S8  G  1  1 

with  positive  and  negative  values  of  7  indicating  countercurrent  and  cocurrent  flows  respectively. 
For  the  case  where  02£2  is  much  smaller  than  4a2?r2  (as  derived  in  eq  18),  the  temperature  distri¬ 
butions  computed  from  eq  18  and  19  should  be  comparable. 


Variable  flow 

The  one-dimensional  pressure  variation  within  the  snow  layer  due  to  atmospheric  changes  of  a 
sinusoidal  nature  can  be  written  as 


where  a  =  kpQ/e 

p0  =  zero  pressure  within  the  snow  layer  at  t  =  0 
e  -  porosity 

k  =  permeability  defined  by  Darcy’s  law;  i.e. 


P'  =  the  difference  between  p2  -  p§ 
P=Po  l1  -M72p0)5in  wf] 

A'/2p0  =  amplitude  of  the  pressure  waves 
u>  =  wave  period. 


(22) 
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Equation  21  is  solved  with  the  following  initial  and  boundary  conditions: 


P'fy.OJ  =0  (23a) 

/5'(0,r)  =0  (23b) 

P\ 8,  t)  -  p 1  -  p\  (23c) 

withp  -  pQ  |1  +  (A'/2p0)  sin  cur)).  Equation  23c  becomes 

P'( 8,  t)  a  A  sin  cat  (24) 


by  neglecting  the  term  A  '2/4  iin 2  uit  and  with  A  =  A  'pQ.  The  solution  of  eq  21  with  the  conditions 
of  eq  23  (Carslaw  and  )aeger  1959)  can  be  expressed  as 


_  fcosh  2 r\y  -  cos  2rjyt 
[cosh  2 pH  -  cos  2 r?8  j 

exp  (-orr2rr2/82) 


sin  (cot  +  ?>)  +  2rro  £ 
n=  1 


r?(-1)n+1to82  nrry 

— — -  sin  -f- 

o2rr27r2+a)2C2  ® 


(25) 


where  r?  -  (co/2a),/j  and  <p  is  defined  as 
0  A  g  sinhr?8(1+/)  ‘ 

Rewriting  eq  25  and  considering  the  steady-state  case,  it  becomes 

p2~po  _  [cosh  2ny- cos  2ryl  . 

4  [cosh  2rj8  -  cos  2t?8  J  Sm  ^  ^  ' 

To  evaluate  the  velocity  through  the  surface,  the  pressure  gradient  is  found  from  eq  27  as 


2p(dJ/ 


:) 


2r?  (sinh  2ry  +  sin  2 qy) 


4  \  9y/  (Cosh  2q8  -  cos  2ni)Vl  (cosh  2rjy  -  cos  2 

.  f  cosh  2»iy  -  cos  2Tjy}  ^  d</> 

+  — rrT - cos  («f  +  0)  "T" 

Icosh  2tj8  -  cos  2tj8  J  dy  ■ 


sin  (o of  +  0) 


(26) 


(27) 


(28) 


The  value  of  dtpjdy  at  y  =  8  is  found  from  eq  26  and 

d<t>  T?(sin  q8  cos  t;8  -  sinh t?8  cosh  pi) 
dy  (sinh  tj8  cos  tj8)2  +  (cosh  tj8  +  sin  r?8)2 

Substituting  eq  29  into  28  and  simplifying,  the  value  of  dp/dy  aty  =  8  is 

dp  _  4q  12( sin  h2qg  +  sin  2q£) 
tfy  y=e  2d  I  (cosh  2tjB  -  cos  2t?8)  sin 

,  (sin  n8  cos  r?8  -  sin  h  r?8  cos  h  n8)  1 

+  — s - '■ - - 1 - ,-L~  cos  cot  . 

(sinh ??8  cos  178)2  +  (cosh  jj8  sin  t?8)2  J 


(29) 


(30) 
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For  small  values  of  rjC,  eq  30  can  be  much  simplified  and  the  velocity  at  the  upper  surface  (i.e. 
y  -  S)  can  be  expressed  as 


y  =  fi 


A  Lk  . 

-SM2?5' 


sin  t of 


2 


COS  GJ/1 


(31) 


Since  the  value  of  rj  is  generally  rather  small  the  second  term  in  the  bracket  can  be  considered 
insignificant  in  comparison  with  the  first  term.  Therefore  the  velocity  at  the  surface  is 


v 


v 


Ak 

/>« 


sin  cjt  . 


(32) 


By  substituting  A  =  A'pQl  and  p  -  p0  1 1  +  (A'l2p0)  sin  wf] ,  eq  32  becomes 


ft  A’ 

f  sin  u>t  1 

"  * 

(1  +  (A’llp0)  sin  cut  1 

Equation  33  indicates  that  the  air  flow  rate  is  sinusoidal  and  varies  in  the  same  manner  as  the  im¬ 
posed  atmospheric  pressure.  The  results  from  eq  1 7  and  1 8  have  shown  (see  next  section)  that 
the  temperature  within  the  snow  layer,  when  there  is  a  constant  flow  of  air,  would  either  be  greater 
or  smaller  than  that  with  heat  transfer  by  conduction  only,  the  difference  depending  on  whether 
the  fluid  flow  is  a  countercurrent  or  cocurrent  in  relation  to  the  direction  of  heat  flow.  If  the 
imposed  flow  is  of  a  periodical  nature,  the  snow  temperature  would  be  expected  to  oscillate  about 
the  neutral  point.  The  repeated  alternate  changes  in  the  direction  of  heat  and  mass  transfer  account 
for  the  almost  constant  snow  densities  observed  near  the  surface  layers. 


COMPARISON  AND  DISCUSSION  OF  RESULTS 
Experimental 

Temperature  profiles  under  pseudo-steady  state  conditions  were  taken  in  connection  with  the 
study  of  the  effective  thermal  conductivity  of  snow  (Yen  1962).  In  that  study,  snow  fractions 
from  sieve  analysis  were  used  with  densities  ranging  from  0.376  to  0.472  g /cm3  and  a  snow  particle 
nominal  diameter  of  0.22  to  0.07  cm.  To  eliminate  or  reduce  the  radial  heat  transfer,  the  snow 
bed  was  contained  in  a  Dewar  flask,  6.5  cm  i.d.  and  15.24  cm  long,  which  was  made  of  a  double- 
walled  glass  tube,  the  annular  space  being  evacuated  and  the  internal  surfaces  plated  with  silver. 
Saturated  air  at  constant  temperature  was  introduced  into  the  top  of  the  vertical  bed  and  passed 
downward,  counter  to  the  upward  flow  of  heat  energy.  The  heat  source  used  to  obtain  a  uniform 
temperature  across  the  bottom  of  the  snow  bed  was  a  copper  cylinder  immersed  in  a  constant 
temperature  bath  below  the  bed.  Two  sets  of  30-guage  copper-constantan  thermocouples  (7  in 
each  set  and  equally  spaced)  were  placed  along  the  center  axis  and  along  the  wall  of  the  bed,  re¬ 
spectively.  Any  difference  between  corresponding  center  and  wall  temperatures  would  indicate 
a  radial  flow  of  heat  through  the  walls  of  the  Dewar  flask. 

The  constant  temperature  bath  had  a  maximum  temperature  variation  of  ±  0.05°C,  and  minor 
fluctuations  were  smoothed  out  by  the  mass  of  the  large  copper  cylinder.  The  temperature  readings 
through  the  snow  bed  were  taken  at  regular  intervals  until  the  pseudo-steady  state  was  arrived  at*. 
When  this  condition  was  reached,  the  air  flow  rate  was  measured.  To  ensure  that  a  steady  state 


*ln  actuality  a  steady  state  will  never  be  achieved  because  of  the  simultaneous  heat  and  mass  transfers.  Due  to 
the  mass  transfer  of  sublimation  and  subsequent  condensation,  the  snow  mass  was  constantly  redistributed  as 
the  experiments  progressed,  thus  effecting  the  conductive  heat  transfer  or  the  snow  bed  temperature  distribution. 
However,  for  the  small  temperature  differences  imposed  between  the  top  and  bottom  of  the  snow  bed,  the  re¬ 
distribution  of  snow  mass  within  the  snow  bed  was  a  rather  slow  process,  so  that  a  pseudo-steady  state  can  be 
considered  to  have  occurred. 
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Figure  1.  Experimental  data  on  dimensionless  temperature 
ip  as  a  function  of  dimensionless  distance  rj. 


had  been  attained,  all  the  read:ngs  were  taken  again  after  a  period  of  30  minutes.  The  time  re¬ 
quired  to  reach  the  steady  state  varied  with  the  rate  of  air  flow  and  was  about  3  to  6  hours  for  a 
mass  air  flow  rate  of  2.84  x  10~4  to  40.17  x  10'4  g /cm2  s.  Figure  1  shows  some  typical  dimension¬ 
less  temperature  distributions  ( <p  -  T-T0/Tt-T0)  within  the  snow  bed  of  density  0.472  g/cm2  s.  It 
shows  a  significant  deviation  from  the  normally  linear  temperature  distribution.  (A  linear  temper¬ 
ature  distribution  will  occur  only  for  the  case  in  which  the  sole  heat  transfer  mechanism  taking  place 
within  the  medium  is  pure  conduction  without  any  involvement  of  mass  transport;  therefore,  it  is 
expected  that  this  non-linear  distribution  will  hold  even  at  no-through-flow  conditions.)  The  de¬ 
viation  apparently  has  been  enhanced  by  simultaneous  mass  transfer  via  the  repeated  process  of 
sublimation  and  condensation. 

In  this  study,  the  heat  is  conducted  upward  and  the  air  flow  is  downward,  and  it  clearly  demon¬ 
strates  that  the  temperature  of  the  snow  is  closer  to  that  of  the  colder  surface  and  is  considerably 
affected  by  the  rate  of  flow.  The  steep  temperature  gradient  in  the  warm  region  in  Figure  1  is  a 
result  of  the  continued  depletion  of  water  vapor  from  a  warm  region  and  deposition  in  a  cold  re¬ 
gion  (prior  to  the  onset  of  the  experiment,  the  snow  bed  had  uniform  density  and  temperature). 

This  metamorphic  phenomenon  and  its  effect  on  the  redistribution  of  the  snow  mass  have  been 
reported  by  Kondrateva  (1945),  Yosida  (1950),  de  Quervain  (1963)  and  Yen  (1963).  If  the  steep 
temperature  gradient  is  sustained  for  a  long  duration,  the  continued  transport  of  water  vapor  will 
eventually  lead  to  the  formation  of  faceted,  cohesion  less  cyrstals  in  the  warmer  area.  A  surface 
layer  of  nearly  constant  snow  density  observed  in  the  Arctic  and  Antarctic  was  duplicated  in  the 
laboratory  experiments  and  is  believed  to  be  caused  by  the  alternate  upward  and  downward 
water  vapor  movement  associated  with  diurnal  temperature  change.  Three  sets  of  data  on  counter¬ 
flow  of  heat  and  fluid  (heat  flow  downwards  and  fluid  flow  upwards,  just  the  reverse  of  the  direc¬ 
tion  in  the  snow  studies  reported  by  Yen  (1962))  in  consolidated  porous  media  (sandstones)  are 
included  in  Figure  1  [Adivarahan  1961 ) .  Adivarahan’s  experiments  were  also  conducted  until 
steady-state  conditions  were  reached.  However,  in  his  study,  no  complications  arose  from  the 
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simultaneous  mass  transfer  via  the  process  of  sublimation-condensation  under  a  constantly  changing 
temperature  gradient.  One  set  of  the  data  was  taken  with  sandstone  135  with  a  porosity  of  0.195, 
a  mean  pore  size  of  9.967  x  10'3  mm,  and  a  nitrogen  flow  rate  of  G  =  73.37  x  10‘4  g/cm2  s.  The 
other  two  sets  were  taken  with  sandstone  of  porosity  0.267,  mean  pore  size  of  1.88  x  10*2 
mm,  and  with  nitrogen  and  carbon  dioxide  flow  rates  of  25.23  x  10'4  and  107  x  10*4  g/cm2  s  re¬ 
spectively. 

Contrary  to  the  case  of  air  flow  through  snow,  the  temperature  distribution  in  the  sandstone 
did  not  deviate  much  from  the  linear  distribution.  In  the  case  of  B0o |d,  it  had  a  greater  porosity, 
and  though  the  nitrogen  flow  rate  was  lower,  it  produced  a  temperature  distribution  almost  identi¬ 
cal  with  sandstone  1 35  with  about  a  flow  rate  three  times  greater.  Since  B0ol(j  is  of  greater  por¬ 
osity,  the  intricate  fluid  velocity  is  lower  in  comparison  with  sandstone  135,  but  the  increase  in 
the  heat  transfer  area  per  unit  volume  of  the  medium  counters  the  effect  of  the  more  intensive 
turbulent  characteristic  in  sandstone  135.  In  the  case  of  B0o |d,  with  carbon  dioxide  flowing  at  a 
rate  about  four  times  greater  than  the  nitrogen,  there  exists  a  much  smaller  difference  (in  compar¬ 
ison  with  the  case  of  air  flow  through  snow)  in  the  dimensionless  temperature  distribution,  indi¬ 
cating  the  more  turbulent  nature  of  flow.  From  these,  as  well  as  the  studies  on  snow,  it  is  clear 
that,  irregardless  of  flow  direction,  the  temperature  is  always  closer  to  that  of  the  colder  surface 
of  the  medium  as  long  as  there  is  counterflow  of  heat  and  fluid. 

Theoretical 

For  the  case  when  02C2/4a2rr2  (a  term  in  the  denominator  of  eq  16)  is  much  less  than  one  and 
for  the  case  of  steady  state,  cq  18  can  be  used  to  compute  the  dimensionless  temperature  distribu¬ 
tion  0.  In  the  case  of  the  snow  experiments,  the  snow  sample  has  a  length  of  1 5.24  cm,  and  the 
value  of  p  is  evaluated  from  eq  9  to  be  /3  =  1.37586.  ae  is  computed  from  eq  8  by  using  Ke  values 
given  by  eq  1 .  For  the  flow  rate  of  G  =  2.84  x  10~4  g/cm2  s,  the  ratio  of  0282  to  4ar2rr2  was  found 
to  be  0.021 ,  a  value  much  less  than  one.  Therefore  for  G  =  2.84  x  10'4  g/cm2  s,  eq  18  becomes 
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Figure  2.  Comparison  of  experimental  data  with  theoretical  predictions. 


Table  1.  Dimensionless  temperature  dis¬ 
tributions  in  a  snow  bed  with  counter- 
current  0  and  cocurrent  flow  of  heat 
and  fluid  (G  =  2.84  x  10*4  g /cm2  s). 


Y/e 

eq  34 

eq  35 

eq  34 

eq  35 

0 

0 

0 

0 

0 

0.1 

0.14S0 

0.1450 

0.0644 

0.0645 

0.2 

0.2776 

0.2775 

0.1349 

0.1352 

0.3 

0.3988 

0.3986 

0.2120 

0.2124 

0.4 

0.5094 

0.5093 

0.2964 

0.2970 

0.S 

0.6106 

0.610S 

0.3889 

0.3895 

0.6 

0.7031 

0.7029 

0.4901 

0.4906 

0.7 

0.7876 

0.7876 

0.6008 

0.6014 

0.8 

0.8648 

0.8648 

0.7220 

0.7225 

0.9 

0.9338 

0.9381 

0.8591 

0.8601 

1.0 

1.000 

1.000 

1.000 

1.000 

For  the  same  flow  rate,  eq  19  reduces  to 

T~T0  =  1  -  exp  (±  Q.059y)  (35) 

0  re-T0  1  -  exp  (±  0.059C)  ' 

As  shown  in  Figure  2  and  Table  1,  for  a  pore  velocity  of  approximately  vt  =  0.48  cm/s.  the  cal¬ 
culated  values  of  0  and  0'  from  eq  34  and  35  are  almost  identical,  indicating  the  validity  of  eq  34 
at  least  up  to  this  pore  velocity  (ve  =  G/epa,  where  e  is  the  porosity  e  =  1  -  pjp\,  and  Ps and  Pi 
are  densities  of  snow  and  ice  respectively).  As  indicated  in  Figure  2,  for  cocurrent  flow  of  heat 
and  fluid,  the  temperature  of  the  medium  is,  on  the  contrary,  closer  to  that  of  the  warmer  end  of 
the  sample.  The  0  values  evaluated  from  the  experimental  results  are  shown  in  the  figure;  a  con¬ 
sistent  variation  from  the  values  computed  from  eq  34  and  35  seems  to  occur  for  G  =  2.84  x  10'4 
cm/s.  However,  the  comparison  between  the  experimental  results  and  computed  values  of  0  from 
eq  35  for  G  =  1 3.1 7  and  40.1 7  x  10'4  g/cm2s,  respectively,  appears  in  good  agreement  for  a  greater 
portion  of  the  sample  lengths.  (Equation  34  cannot  be  used  for  comparison  because  in  these 
cases  the  values  of  02$2  and  4a2 n2  are  of  the  same  order  of  magnitude.) 


CONCLUSIONS 


From  this  investigation,  it  can  be  concluded  that  for  counterflow  systems  and  for  those  involving 
vapor  diffusion  through  the  process  of  sublimation  and  condensation,  the  dimensionless  tempera¬ 
ture  0  is  always  greater  than  the  value  from  the  linear  temperature  distribution  for  no-through- 
flow  and  no-moisture-transport  situations.  The  temperature  gradient  is  found  to  be  much  steeper 
at  the  warm  region  than  at  the  colder  regions  and  the  accompanying  moisture  transfer  processes 
under  a  temperature  differential  significantly  influence  the  redistribution  of  snow  mass,  i.e.  a  net 
loss  in  the  warm  zone  and  a  net  gain  in  the  cold  zone.  The  extent  of  this  mass  redistribution  de¬ 
pends  on  the  overall  temperature  difference  imposed,  the  temperature  level,  and  the  mass  flow  rate. 

For  counterflow  systems  without  accompanying  simultaneous  vapor  transport,  the  nonlinearity 
in  the  dimensionless  temperature  0  is  remarkably  reduced.  For  a  cocurrent  flow  system,  the  re¬ 
verse  of  the  observed  temperature  distribution  exhibited  in  a  counterflow  system  is  true;  i.e.  tem¬ 
perature  gradients  are  steeper  in  the  colder  zones,  and  temperature  gradients  increase  very  slightly 
in  the  warm  zones,  for  a  great  portion  of  the  sample  depth. 
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